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I. INTRODUCTION 

Recently multi-dimensional quantum systems (quan- 
tum dits, or qudits) attract much attention in context 
of quantum information and communication. It is partly 
caused by fundamental aspects of quantum theory since 
usage of qudits allows one to violate Bell-type inequalities 
strongly than with two dimensional systems (qubits) pj . 
Much interest in qudits also comes from the application 
point, especially from applied quantum key distribution. 
Multilevel systems are proved to be more robust against 
noise in the transmission channel, although measurement 
and preparation procedures of such systems seem to be 
much more technically complicated than in the case of 
qubits. Different aspects of the security of qudits-based 
protocols have been analysed including those related 
to the reduced set of bases • Lately a proof-of- principal 
realization of a QKD protocol with qudits 0] and with en- 
tangled qutrits (D = 3, where D is dimensionality of the 
Hilbert space) [fj have been demonstrated. For the last 
several years elegant experiments were performed where 
different kinds of optical qudits were introduced. Most 
of them are based on states of light emitted via spon- 
taneous parametric down-conversion (SPDC). There are 
qudits obtained using spatial degrees of freedom: with 
either single photons Q or twin photons [||; orbital an- 
gular momentum of photons 0; time-bin tecnique 
multi-armed interferometers |9j; postselection of polar- 
ization qutrits from four-photon states [id. and polar- 
ization states of single- mode biphotons [TlJ . 

It seems that among manifold manners of qudits prepa- 
ration only the method based on polarization states of 
single-mode biphotons admits preparation of any qutrit 
state (pure or mixed) on demand within one set-up and 
guarantees complete control over the state with high ac- 



curacy |llill2j . However, this method does not allow to 
create entangled qutrits. Another limitation of polariza- 
tion qutrits relates to the fundamental lock for their basic 
states transformations using SU(2) optical elements, like 
retardation plates, rotators etc. Meanwhile it is these 
transformations that would be extremely useful for quan- 
tum communication purposes. 

In this paper we present the results of experimental 
preparation and measurement carried out with polariza- 
tion based ququarts or quantum systems with dimen- 
sionality D = 4. The paper is organized as follows. In 
Sec. II we discuss the main properties of ququarts based 
on the polarization state of two-photon (biphoton) field. 
Such operational notions as coherence matrix, polariza- 
tion degree are introduced. The criterium of separability 
for qubits forming ququarts is discussed as well. Sec. Ill 
is devoted to different experimental implementations of 
biphoton-ququarts preparation and their measurement. 
In Sec. IV we consider a specific set of ququarts which 
can be used in practice for quantum key distribution. 



II. POLARIZATION QUQUARTS AND THEIR 
PROPERTIES 

A. Polarization states of biphoton field 

Biphoton field generated via spontaneous paramet- 
ric down-conversion process is represented by a coher- 
ent mixture of two-photon Fock states and the vacuum 
state H3|: 



* = \vac H^E^I 1 ^)' (!) 
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where |l,-,l,-») denotes the state with one (idler) pho- 
ton in the mode ki and one (signal) photon in the mode 
k s . Since squared modulus of the coefficient Fg g gives 

probability to register two photons in modes ki and k s it 
is called the biphoton amplitude 0|. The pure polar- 
ization state of down-converted light field, which is often 
refered as two-photon field (or biphotons), can be written 
as follows: 
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amplitudes. The first (n) and second (to) place in kets 
corresponds to number of distinct photons with definite 
(either horizontal H or vertical V) polarization, with to- 
tal photon number m + n = 2. If down converted photons 
have equal frequency and momentum, then a ququart 
state (0 converts to a qutrit state, i.e. middle terms 
in (|5J| become indistinguishable: k s « ki, ui s w uji and 
lo s + uii = lo p , where uj p is the laser pump frequency. In 
this case any arbitrary pure polarization state of bipho- 
ton field can be expressed in terms of three complex am- 
plitudes Ci , c' 2 , and c' 3 : 



|c')=ci|2,0)+4|l,l)+4|0,2), (3) 



with c- = |c£| exp{i(^ 4 }, J2 \ c i\ 2 = 1- The vector |c') = 

(c? ly c* 2 , c' 3 ) represents a three-state system or a qutrit. 
Figure 1 presents the photograph picture of two dimen- 
sional SPDC spectrum taken in coordinates wavelength- 
angle. By convention the lcft-and-uppcr side in re- 
spect with doubled pump wavelength belongs to "sig- 
nal" range whereas the right-and-down side belongs to 
"idler" one. The dashed lines limit the central part of 
the spectrum which corresponds to the frequency de- 
generate and collinear regime when photons forming the 
biphoton have approximately the same wavelengths and 
propagate along the pump. To select this part pinholes 
and/or narrow-band filters are typically used. It is that 
part of biphoton spectrum, that contributes to qutrits. 
In order to distinguish between middle terms in @ one 
must induce the distinguishability between down con- 
verted photons cither in frequency, momentum or detec- 
tion time. In experiments described in this paper we 
chose collinear non-degenerate regime of SPDC so twin 
photons that form a biphoton were having different fre- 
quencies and propagate along the same direction. The 
appropriate two-dimensional spectrum is shown in Figure 
2. To take this spectrum it is enough to tilt the crystal by 
a small angle with respect to the orientation that is used 
for degenerate regime. Two dashed rectangles indicate 
the angular-frequency ranges of spectrum contributing 
to ququarts. The state can be thus rewritten in the 
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FIG. 1: Photograph of lithium iodate (LHOs) frequency- 
angular spectrum taken with helium cadmium laser operated 
at 325nm. Angle between crystal optical axis and laser pump 
is 59.22°. Dashed lines select the part of the spectrum (fre- 
quency degenerate, collinear regime of SPDC) contributed to 
biphotons-qutrits. Vertical lines at the bottom belong to the 
spectrum of neon which serves as a frequency reference. 
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FIG. 2: Non-degenerate regime of SPDC. Angle between crys- 
tal optical axis and laser pump is 58.97°. Dashed lines select 
the part of the spectrum contributed to ququarts. 



form: 



c 1 \H 1 ,H 2 )+c 2 \H 1 ,V 2 )+c 3 \V 1 ,H 2 )+c i \V ll V 2 ) 



Here \H l{2) ) = a\ 



'1(2)' ^1(2) are t ne creation operators of down-converted 
photons with central wavelength Ai(A2) in horizontal and 
vertical polarization mode. The alternative way to repre- 
sent a ququart might be by introducing distinguishabil- 
ity between spatial modes while keeping the same central 
wavelengths of photons: k s ^ ki, lu s w u>i. The detailed 
description of polarization properties of such a state of 
two qubits has been done by James and co-authors [Tlj . 
The case of ququart with two frequency modes attracts 
much attention because for some class of the tasks it is 
convenient to operate with states in single spatial mode. 



(2) 



|Vl(2 



(4) 

= b^ 2 j\vac), where 



3 



B. Is ququart a separable or entangled state of two 
qubits? 

Generally, the state J3J can not be written as a direct 
product of polarization states of two photons: 



l*i> 



{pi\H 1 )+q 1 \V 1 ))®{ P 2\H 2 )+q 2 \V 2 )) (5) 



That is why it is useful to write down the condition when 
@ becomes a separable state. The reduced density ma- 
trix of the subsystem (second photon) can be found by 
tracing the joint state over all states of the first photon: 



P2 = s Pl mm 



Cl 



\C3\ CiC 2 +C 3 C 4 



C2C!+C4C3 \C 2 \ 



C4 



(6) 



It is straightforward to show that the eigenvalues of © 
take the meanings Xi f2 = 0, 1 and entropy of each sub- 
system drops to zero iff 



Cl C4C4 -I-C3C3 C2 C2 — C2 C3 C4 — C4 C3 Cl C2 = |ciC4 — C2C3I 
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(7) 

Thus if C1C4 = C2C3 holds, the state (@J reduces to factor- 
ized one with definite states of each qubit l^i), l 1 !^) and 
vice versa. In other words J7J) plays the role of criterium 
for biphoton-based ququart to be a separable state of a 
couple of polarization qubits belonging to two different 
modes. Sub-indexes in Q designate either frequency or 
spatial modes. 

As examples of different ququarts we refer to the prod- 
uct states that might be generated from a single non- 
linear crystal via SPDC: 



typel: \HiH 2 ) 
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and maximally entangled Bell-states: 
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At the end of this section we would like to mention 
that biphoton-ququart is not supposed to be mapped 
on Poincare sphere as a pair of points like it occurs for 
biphoton-qutrits [lit Il7j . Indeed it is not because the 
state (0} can not be factorized in the general case. How- 
ever, if (JJJ holds the sub-set of separable qubit states 
admits its visual representation on Poincare sphere. 



C. Second order of the field: Stokes parameters 



of Stokes operators averaged over the state Although 
the description of light polarization can be introduced 
only for the quasimonochromatic plan e waves, it is possi- 
ble to use P-quasispin formalizm [l8j to describe the po- 
larization of arbitrary quantum beams with n modes, fre- 
quency or spatial. For two-frequency and single-spatial 
mode field, the Stokes parameters will contain time and 
frequency dependent terms exp(i(ujj —u>f-)t) that describe 
" beats" of frequency modes and have no connection with 
light polarization. However, these terms vanish if one 
considers the finite detection time that allows to clas- 
sically average these "beatings". So in the case of two 
frequencies, the Stokes operators are given by the sum of 
corresponding operators in two modes 
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2i?e(c*(c 2 +c 3 ) +c 4 (c* +c*)); 
(S3) = (a\bi + a\b 2 - b\ai - b\a 2 ) = 

(4 1} > + (sP) = 

2Im(c\(c 2 + c 3 ) + c 4 (4 + C3O). 

The same definition applies to the second wavelength 
A2 and we take into account that these operators do not 
commute for different frequency modes. 



D. Fourth order of the field: Coherence matrix 

Polarization properties of biphoton-ququarts are de- 
termined completely by the analogue of the coherence 
matrix. It is a matrix consisting of ten fourth-order mo- 
ments of the electromagnetic field. An ordered set of 
such moments can be obtained using the direct product 
of 2 x 2-coherence matrixes for both polarization qubits 
forming biphoton: 



K 4 = (K 2 )i <g> (K a )i 



A E F G 

E* B I K 

F* I* C L 

G* K* L* D 



(11) 



where (^2)1,2 are 2 x 2-coherence matrixes for single 
photons with different frequencies marked with indexes 
J = 1,2: 



(K 2 )j = 



(bfa) 



(12) 



Polarization properties of a pure state (@J can be de- 
scribed by Stokes parameters which are the mean values 



The diagonal elements of (|llfl are formed by real mo- 
ments, which characterize the intensity correlation in two 
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polarization modes H and V: 



A = (a\a\a\a 2 ) 
C = (b\alb ia2 ) 



| Cl | 2 , B={a\b\ ai b 2 ] 
|c 3 | 2 , D = {b\b\bxb2) 



|c 2 | 



Nondiagonal moments are complex: 



(13) 
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= {a\b\bib 2 ) = c|c 4 , 
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C2C3, 

C 3 C 4- 



(14) 

The averaging in (|13l I14|l is taken over the state Q). 
The polarization (reduced) density matrix of ququart co- 
incides with coherency matrix K4. 



E. Polarization degree and transformations of 
ququarts 

The polarization degree is given by 



P 



1 E (s k ) 
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This definition of polarization degree is just generaliza- 
tion of the commonly used classical one. It differs from 
the definition suggested in 0] , where it serves as a wit- 
ness of the state pur ity. In the case of polarization-based 
qutrit states [Til Il2| , the polarization degree is 
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with c'j = C\ 1 \[2c! 2 = C2 = C3,Cg = C4. This quan- 
tity was invariant under unitary polarization transforma- 
tions, which is reasonably straightforward [2(j. Indeed it 
was impossible to prepare all demanded pure states that 
can be used, for example, in QKD protocols, unless one 
uses interferometric schemes with several nonlinear crys- 
tals [Til ], or introduces losses. In particular, there is no 
way to transform the basic qutrit state = \V, V) 

with P = 1 into orthogonal state \^' 2 ) = \H,V) with P 
= using retardation plates [2l| . However, in the case of 
polarization ququarts, this quantity is no longer invari- 
ant. Indeed the polarization degree l|15|) can be rewritten 
as 



(i)e(2)\ 
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(17) 

In (I17|) we used trivial links between coherence matrices 
and Stokes parameters |23[ . 



The expression in the round brackets is invariant under 
unitary transformations as well as the denominator. The 
second term under the square root in numerator can be 
expanded in terms of fourth-order moments: 



{1)Q i 2) ) = 



{TrK 4 -2(ReI + B + C)}. (18) 
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Since unitary transformations keep number of photons, 
then TrKu is invariant. At the same time it is easy to 
prove that combination of moments Rel + B + C is not 
an invariant, so the whole expression (|17f) changes under 
SU(2) transformations. As a consequence, the polariza- 
tion degree changes by applying local unitary transfor- 
mations in each frequency mode. This can be achieved 
by using dichroic polarization transformers, which act 
separately on the photons with different frequencies. For 
example, to transform the state 



1*4 



\V Xll V X2 ) => |*2) = \H Xl ,V X2 



(19) 



one needs to use the retardation plate which serves as a 
half wave plate at Ai and as a wave plate at A2-For the 
general case when ququart is not a product state of two 
qubits the transformation matrix G: 



G\V) 



(20) 



can be found simply by making use of the Heizcnbcrg 
picture. We use the equivalent representation of (0J using 
creation and annihilation operators: 



\Ht) m = [ Cl a\al + c 2 a[bl + c 3 b[al + c 4 blbl}\vac) . (21) 

SU(2) transformation between input aj, bj and output 
a'j, b'j operators is given by: 



where 



a'j = tjdj + rjbj 
b'j = -r*aj+t*bj, 



cos 5j + i sin Sj cos 2a 
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Here tj and Tj are the amplitude transmission and reflec- 
tion coefficients of the waveplate, Sj is its optical thick- 
ness, h is the geometrical thickness, aj is the orientation 
angle between the optical axis of the plate and vertical 
direction. Lower indices j = 1 , 2 numerate the frequency 
modes of photons. Substituting l|22|) into l|21|) one can im- 
mediately find that the unitary transformation on state 
(0} is given by 4 x 4 matrix which is obtained by a direct 
product of two 2x2 matrices describing the transforma- 
tion performed on each photon pol |: 



G : 
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It is important to note that the optical thickness 5 
depends on the wavelength, so parameters t, r which de- 
termine the result of transformation differ for photons 
forming biphoton-ququart. Thus for the ququart trans- 
formation (without taking non-essential phase term into 
account) considered above Ijl9(l the matrix has the form 



10 
1 
10 
10 



(25) 



The same dichroic plate performs unitary transforma- 
tions between polarization Bell states: g]^) = | x f ,± - ) ), 
and g|4' ± ) = | < I )± ), where and l 1 ^) are represented 
by ququarts with C\ = ±C4 = and C2 = ±03 = -7= 
correspondingly. Similar transformations with frequency 
non-degenerate biphotons prop agating in single spatial 
mode have been realized in [22 ■ 



III. EXPERIMENTAL IMPLEMENTATION 

In this section wc consider a several methods of 
ququart state preparation and their measurements. 



A. Preparation 

In general to prepare arbitrary ququart state (0J it is 
necessary to use four nonlinear crystals arranged in such 
a way that each crystal emits coherently one basic state in 
the same direction. But in particular cases reduced set 
of crystals is quite enough to generate specific ququart 
states which can be used in practice. In experiments de- 
scribed below we used the following methods to generate 
sub-set of ququarts. 



2. Entangled states 

To prepare the ququarts with C1C4 7^ C2C3 it was suf- 
ficient to use two crystals like it was done when fre- 
quency non-degenerated Bell-states have been generated 
M (Fig.4). 

In order to introduce a phase shift 9314 between hori- 
zontally and vertically polarized type-I biphotons two 1 
mm quartz plates (QP1,2) that can be rotated along the 
optical axis were used. By tilting the plates the family 
of ququarts was generated: 



1*1/(0)) 



|ci| 





(27) 



c 4 e 



where real amplitudes | ci | , | C4 1 were controlled by the 
half wave plate inserted into the linearly polarized pump 
beam. Some remarkable interference phenomena with 
|\I/j/(0))-like states have been observed. For example, 
when the phase of Ai-photon is varied, then modulation 
with the same wavelength is observed in coincidence rate, 
while changing the phase of A2-photon it is this modu- 
lation that is manifested in coincidences [24|. Another 
effect related to biphotons-ququarts that is widely dis- 
cussed in quantum optics is "hidden polarization". Po- 
larization transformations of frequency non-de gen erate 
Bell-states experimentally have been studied in [22| . 



pump 




pump 



1. Separable states 

If nonlinear crystal generates one of the basic state 
© then applying the transformation (|24f) one can get 
the sub-set of ququarts which are the product states of 
pair polarization qubits (Fig. 3). We used the basis state 
IV1V2), so the states generated in this way have the fol- 
lowing structure: 



!*/(")> 




(26) 



Since coefficients rj and tj depends on the the orientation 
angle aj there is a simple way to transform the state by 
rotating the retardant plate. As we will show below, a 
usage of a single crystal allows one to prepare the whole 
sub-set of ququarts which can be used for quantum key 
distribution. 



FIG. 3: Preparation of states VP/. The nonlinear crystal gen- 
erates the basis ququart state(s) via type-I SPDC. Then the 
state is transformed by the wave plate QP whose parameters 
depend on its orientation angle a. 




BBC.) 



QIT QP2 



FIG. 4: Preparation of states ^n. Two nonlinear crystals 
generates pair of the basis ququart states via type-I SPDC. 
The relative phase between states ipu is controled by two 
quartz plates QP1, QP2. The state amplitudes are con- 
trolled with the help of half-lambda plate. 
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B. Measurement 

Basically it is necessary to perform D 2 independent 
measurements for complete reconstruction of the arbi- 
trary qudit density matrix. So for D = 4 the total num- 
ber of measurements is equal to sixteen. Since the only 
way to measure biphoton is to register a coincidence of 
photocounts, we chose a Brown- Twiss scheme supplied 
with polarization elements to perform projective mea- 
surements. A coincidence click that occurs when signals 
from two detectors fall into the coincidence window is 
considered a registered event. Since the coincidence click 
appears at the output with some probability, a statisti- 
cal treatment of the outcomes becomes very important. 
Another point that we would like to notice is that ac- 
cording to Bohr's complementarity principle, there is no 
way to measure all moments <|13l H4|) at the same time, 
dealing with a single ququart only. So to perform a com- 
plete set of measurements one needs to generate a lot of 
representatives of a quantum ensemble. 

In order to measure the ququarts we developed two 
protocols. 



NPBS IF /A fj\ V 




FIG. 5: Measurement block for Protocol 1. The Brown- 
Twiss scheme for measuring intensity correlation between 
two polarization modes. After spatial separation at the non- 
polarizing beam splitter (NPBS), signal (A2) and idler (Ai) 
photons propagate through the quarter- (QWP1,2) and half 
(HWP1,2) wave plates , polarizing prisms (PP1,2) in two 
channels. The interference filter (IF) placed in the upper 
arm and centered at Ai selects frequencies modes of ququart. 
Finally, photons are registered by detectors (Dl,2). The co- 
incidence rate from the output of the coincidence circuit (CC) 
is proportional to the fourth moment in the field (R). 



1. Protocol 1. 

The idea of the first method is to split the ququart state 
\^in) into two frequency modes and to perform transfor- 
mations independently on each of the photon from a pair 
(Fig. 5). Usage of a dichroic beamsplitter allows one 
to separate frequency modes in space and reduces this 
protocol to that one developed in for two spatially 
separated polarization qubits. We called this method 
as "frequency selective" because frequency mode sepa- 
ration is basically applied to ququart before subjecting 
it to polarization projective measurements. The polar- 
ization transformations can be done using polarization 
filters placed in front of detectors. Each filter consists 
of a sequence of quarter- and half-wave plates and a po- 
larization prism, which picks out vertical polarization. 
In our experiments wc kept the protocol developed in 
|15| unchanged, although there are many other configu- 
rations of polarization filters leading to the matrix K4 
reconstruction. Instead of dichroic beamsplitter we used 
a (50/50)% polarization insensitive beamsplitter and a 
narrowband filter centered at Ai was put in front of the 
first detector. In such a way we postselected the desir- 
able events in registering coincidences. Approximately in 
quarter of trials, the photons with central wavelength Ai 
forming a biphoton are going to the first detector, while 
the conjugated ones (with central wavelength A2) are go- 
ing to the second detector. In the remaining cases, the 
events are not registered because they either do not con- 
tribute to coincidences or corresponding photons do not 
pass through the interference filter. 

In the Heisenbcrg representation the polarization 
transformation for each dichroic beam-splitter output 



port is given by: 








5)1^/2(^=^/2,0) 




(28) 



We designate with j = 1,2 different output ports of 
beamsplitter which correspond to different wavelengths 
Ai ; 2- Four 2x2 matrixes in the right-hand side of 128|) 
describe the action of the non-polarizing beam-splitter, 
A/4-, A/2- plates and vertical polarization prism on the 
state vector of each photon; 



D 



Aj/2,A J /4 



where rj and tj - are the coefficients introduced in Eq. 
(23), so for a Aj/4-plate {Sj = n/A), 



1 i 

= -j=(l+icos2xj),r Xj /4 = -^sm2 Xj (29a) 



and for a Aj/2-platc (Sj = 7r/2) 



t\ 3 /2 =ico8(26 j ),r Xj /2 = isin(20,-). 



(29b) 



Thus, there are four real parameters (two for each 
channel) that determine polarization transformations. 
Namely, these parameters are orientation angles for two 
pairs of wave plates: 0i,Xi,02)X2- Also wc would like 
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to notice that there is another parameter that affects 
strongly on result of polarization transformations. It is 
a wavelength(s) of down-converted photon(s). For exam- 
ple varying the frequency spectrum of SPDC by tilting 
a crystal which generates biphotons one can significantly 
change the transformed state 1|20|) using fixed parame- 
ters of the wave plates such as geometrical thickness h 
and its orientation a. This property makes biphoton- 
based ququarts much more flexible to be controlled than 
biphoton-qutrits and allows one to choose convenient 
regimes for operations with them. 

As it was mentioned above, the output of the Brown- 
Twiss scheme is the coincidence rate of the pulses coming 
from two detectors D\ and D 2 . The corresponding mo- 
ment of the fourth order in the field has the following 
structure: 



R h2 cx (fc'iW&'i&'a) = R(e 1 , X i,02,X2) 



(30) 



In the most general case this moment contains a linear 
combination of ten moments (|13l I14f> forming the matrix 
So the main purpose of the state reconstruction pro- 
cedure is extracting these moments from the set-up out- 
comes by varying the four parameters of the polarization 
Brown-Twiss scheme. Consider some special examples, 
which give the corresponding lines in the complete pro- 
tocol introduced below (Table I). The measurement of 
first four moments (|13f) is trivial procedure. For instance 
the third line in the (Table I) corresponds to selection 
the basis state IU1V2). The other three upper lines corre- 
spond to the measurement of other basis states \HiH 2 ), 
\H\V 2 ) : and \V\H 2 ). Remaining lines of protocol show 
how to measure one of the complex moments 1)14)1, To 
measure the real part of the moment E, let us set the 
wave-plates in the Brown-Twiss scheme in the following 
way. 

The first channel: 



A/4: xi =0°,£> 



A/4 



1 

71 



l + i 






1 - i 



X/2:6 1 =45°,D x/2 

The second channel: 

A/4: X 2 = -45 ,£ A/4 



i 

1 



J_ / 1 -i 



(31a) 



(31b) 



(32a) 



annihilation operators it is easy to get the final moment 
to be measured: 



R = (^16^6^1 6 2 |*) = 1/8(4 + C + 2Re£T). 

A complete set of the measurements called the tomog- 
raphy protocol is presented in Table I. Each row corre- 
sponds to the setting of the plates to measure the mo- 
ment highlighted in the sixth column. The last column 
corresponds to amplitude of the process (see below) . 

Amplitude of the quantum process, corresponding to 
passing down-converted photons through the measure- 
ment set-up for each configuration of wave plates is de- 
scribed by the following equation: 

AI l/ (5 1 ,8 k ,5 2 , (pi) = ^[aia 2 ci + a 1 b 2 c 2 + b x a 2 c z + b x b 2 c^. 



where 



ai = 
a-2 = 
bi = 
b 2 = 



-rx/ 2 {0i)t* x/i {xi) - t\/2(&i)r x /i{xi), 
-r\/2(0 2 )tl/4(x2) - *A/2(02)r- A / 4 (X2), 
-rx/ 2 (0 1 )r* x/i (xi) + tx /2 {6i)t x/i (xi), 
-r\/ 2 (0 2 )r* x/4 (x2) + *A/2( 2)*A/4(X2)- 



(33) 



(34) 



The complete reconstruction of the input state |<F m ) 
was performed according to the algorithm considered in 
our previous work [T^ . 

Concluding this section we would like to stress that 
if usual frequency insensitive (broadband) beamsplitter 
is used and no interference filters are inserted either in 
front of or behind this beamsplitter, then "non-selective" 
method of ququart measurement is possible. In this case 
one does not need to use the interference filter to select 
the wavelength in each channel of Brown-Twiss scheme. 
In other words each detector is allowed to register pho- 
tons with different frequencies and the fourth moment in 
the field to be measured becomes 



Rl,2 OC <(&'! + b' 2 )\ 



Detector2 



(b\ + If 2) 



Detectorl / 



(35) 



instead of (|30|) . In non-selective method the expression 
connecting observable value R± t2 with components of co- 
herence matrix l|13l I14[l is bulky and it is not reasonable 
to utilize it for the state reconstruction. However in par- 
ticular case we applied non-selective method to ququart 
reconstruction procedure and developed the special pro- 
tocol, which is introduced in the next section. 



2. Protocol 2. 



A/2 : 9 2 = -22.5°, D x/2 = -L ( j. J. j . (32b) 

Substituting these matrices into Eq. (28) and taking 
into account the commutation rules for the creation and 



In the second protocol the whole ququart state J2J is 
subjected to the linear transformations using a set of 
the retardant plates |25j |. Then projective measurements 
were applied to the state |\I> m ) after it being transformed 
by the plates. Using the retardant plates with fixed op- 
tical thickness, one can reconstruct the state, varying 
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TABLE I: Protocol 1. Each line contains orientation of the half (6 3l i) and quarter (Xs,i) wave plates in the measurement block. 
Last two columns show the corresponding moment R v and the process amplitude M v {y — 1, ..16). Sign "-" in orientation of 
plates for reflection channel is introduced to compensate 7r-phase shift occurring by reflectance from the beamsplitter. 





Set-up parameters 


Moments to be measured 


Amplitude of the process 


V 


Xs 


0s 


Xi 




R* = \M V \" 


M v 


1. 





45° 





-45° 


A/ 4 


a/2 


2. 





45° 








B/4 


c 2 /2 


3. 














D/4 


c 4 /2 


4. 











-45° 


G/4 


C3/2 


5. 





22.5° 





—45° 


l/8(A + G + 2ImL) 




6. 





22.5° 








l/8(B + D + 21mK) 


^(c 2 -iC4) 


7. 


45° 


22.5° 








l/8(B + D - 2ReK) 




8. 


45° 


22.5° 





-45° 


1/8(A + C - 2ReF) 


275 ( Cl " C3 ) 


9. 


45° 


22.5° 





-22.5° 


1/16(A + B + C + D)- l/8(Im£ + ReL - ImG + Iml + RcK + ImL) 


j[{ci - C 3 ) + j(c 2 - C 4 )] 


10. 


45° 


22.5° 


-45° 


-22.5° 


1/16(A + B + C + D)- l/8(ReF - ReL + ReG + Rel + ReK - ReL) 


3 [Cl + C'2 - C3 - C 4 ] 


11. 





22.5° 


-45° 


-22.5° 


1/W(A + B + C + D) + l/8(ImF + ReL + ImG + Imi + ReL + ImK) 


j[(ci + c 2 ) - i(e 3 + c 4 )] 


12. 





45° 


-45° 


-22.5° 


1/8(A + B + 2ReE) 




13. 








-45° 


-22.5° 


l/8(G + D + 2ReL) 




14. 








-90° 


-22.5° 


l/8(G + L> + 2ImL) 




15. 





45° 


-90° 


-22.5° 


1/8{A + B + 2ImE) 


272^ ~ iC *) 


16. 





22.5° 


-90° 


-22.5° 


1/16(A + B + C + D) + l/8(ImF + ImE - ReG + ReJ + ImL + ImL") 


j[(ci - c 4 ) - i(c 2 + ca)] 




QP1 QP2 



V U^l| PP2 




c c. 



FIG. 6: Measurement block for Protocol 2. After the input 
state is subjected to transformations by wave plates (QP1, 
QP2), the projective measurements selecting vertical polar- 
izations of ququart components were performed. 

the orientation of the plates (Fig. 6). Mathematically 
one needs to derive an equations set connecting real and 
imaginary parts of moments (|13l I14|) with coincidence 
rate. It turns out that the complete set of moments (|13l 
I14f> can be extracted from the projective measurements 
if two different plates were used. 

\* out )kt = G(5 1 ,6 k )G(5 2 ,<pt)\** n ), (36) 

where 0^ and ipi are the orientation angles of the first 

and second plates. The parameters of the plates, i.e. op- 

(1 2) 

tical thicknesses for different wavelengths <5L' and their 
orientations are supposed to be known. The projective 
measurements were realized by means of pair of polariza- 
tion prisms, transmitting vertical polarization and settled 
in front of each single photon detector. So the number of 



coincidences is 

Ru « \{ViV 2 \^out)ki\ 2 - (37) 

In order to reconstruct the state it is sufficient to per- 
form at least four measurements corresponding to differ- 
ent orientations of the second plate (pi and repeat this 
procedure at least four times varying orientation of the 
first plate 9k- But in experiments we performed redun- 
dant number of measurements to increase the accuracy 
and finally used 36 orientation of the second plate. Thus 
total number of measurements /i in this protocol was 
equal to 4 x 36 = 148. Therefore the Protocol 2 includes 
148 lines instead of 16 lines used in Protocol 1. 



C. Experimental setup. 

For generation of biphoton-based ququarts we used 
lithium-iodate 15 mm crystal pumped with 5 mW cw- 
helium-cadmium laser operating at 325 nm. The orien- 
tation of the crystal is chosen at 58.98° with respect to 
the pump wave vector that the down-converted photons 
at Ai = 702nm and A2 = 605nm had been generated. 
For some particular cases we selected biphoton-ququarts 
at wavelengths A 4 = 667nm and A2 = 635nm. Thus 
the either IV702V605) or IV667V635) states were used as an 
initial states. Then, this state was subjected to trans- 
formations done by dichroic waveplates to prepare the 
sub-set of ququarts with C1C4 = C2C3 l|26|) . That sub-set 
of states was used to be reconstructed. In particular we 
used the 0.988-mm (or 0.315-mm ) length quartz plate 
and changed it orientation a, which served as a param- 
eter. Since the thickness of the plate, quartz dispersion 
and orientation a are supposed to be known, we were able 
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to calculate the result of the state transformation with 
high accuracy. When Protocol 2 was applied for ququart 
reconstraction we used the quartz plates with thicknesses 
0.821mm (QP1) and 0.715mm (QP2). Both protocols 
were applied to reconstruct the initial state |^j(a)). In 
the case when sub-set c\ca ^ C2C3 had been generated we 
used two consecutive 1.8 mm thick type-I BBO (beta- 
barium borate) crystals whose optical axis are oriented 
perpendicularly with respect to each other at 36.33° with 
respect to the pump wave vector. An interference fil- 
ter centered either at 702nm or at 635nm and with a 
FWHM bandwidth of 3nm was placed in transmitted 
arm to make a spectral selection of one photon of a pair, 
while the photon with conjugated frequency was sent 
to a reflected arm. As it was already mentioned above 
this scheme is equivalent to that one used in ^5|. The 
only difference between our scheme and developed in [l5j 
resulted from operating with frequency non-degenerate 
collinear regime of SPDC instead of non-collincar, fre- 
quency degenerate regime used in 0. Without loss of 
generality, this scheme also allows to reconstruct any ar- 
bitrary polarization ququart state by registering coinci- 
dence rate for different projections that are done by the 
polarization filters located in each arm. Each filter con- 
sists of a zero-order quarter- and half-waveplate and a 
fixed analyzer. Two Si-APD's linked to a coincidence 
scheme with 1.5 nsec time window were used as single 
photon detectors. 

D. Results and discussion 

We applied both protocols to measure the set of 
states ^>i, generated with \Ve67nmVe35 nm ) transformed 
by quartz plate QP, like it is shown in Fig. 3. We have 
tested the states which were generated for two orientation 
angles of plate QP a = 0°,90°. According to the Proto- 
col 2 four sets of measurements was performed for each 
input state, so totally we performed 148 measurements of 
coincidence rate as a function of orientation angles ipi of 
QP1 and 9k of QP2. The dependence of coincidence rate 
on the orientation <p t for 9 k = 90° (a), 105° (6), 120° (c), 
and 135° (d) is shown in Fig. 7 for the input state corre- 
sponding to orientation a = 20° of the plate forming the 
state. 

The solid curve gives the theoretical behavior of the 
normalized fourth moment, the dots with bars show the 
experimental data. The tables II-IV show the results 
of statistical reconstruction of ququarts \^i) performed 
with root estimation method. The tables show theo- 
retical and experimental four-dimensional state-vectors 
Ctheory, c exp as well as fidelity F defined by 

F = KctheorylCexp)] 2 ■ (38) 

In table II we collected the data measured by Proto- 
col 1 (ququarts with Ai = 702nm and A 2 = 605n,m, 
/i(QP)=0. 315mm). Table III corresponds to measure- 
ments performed by Protocol 2 (ququarts with Ai = 




9(deg.) 



FIG. 7: Dependence of coincidence rate on tilting angle 9 



667nm and A 2 = 635nm, /i(QP)=0.988mm). And fi- 
nally we demonstrate the validity of preparation, trans- 
formation, and measurement procedures performed over 
biphotons-ququarts by reconstruction of the Bell states 
|$ ± ) (Ai = 702?im and A 2 = 605rnn), measured along 
with Protocol 1 and shown in Table IV. As it is seen from 
the tables the accuracy of state reconstruction is a little 
bit lower when Protocol 1 has been applied. It results 
from the fact that Protocol 1 exploits minimal number of 
measurements (16) which is needed for four-state system 
reconstruction. At the same time since redundant num- 
ber of measurements (148) was involved when applying 
the second protocol, the highest accuracy was achieved 
like in the case of biphotons-qutrits ^3 ■ The respectively 
low fidelity for Bell states reconstruction (Table IV) is 
explained by experimental imperfections at the prepara- 
tion stage. Preliminary estimation of the prepared state 
quality can be extracted from visibility displayed during 
interference experiments with these states. For exam- 
ple, typical visibility meaning revealed when space-time 
interference observed was 0.9 [2fi(|. Another factor limit- 
ing the accuracy of the state reconstruction is the finite 
number of events to be registered. In the case of |^ r j)-set 
of states we collected about several thousands of coin- 
cidences during accumulation time whereas in the case 
of |\&//)-set only a few hundreds coincidences were regis- 
tered in total. 



IV. POLARIZATION QUQUARTS IN QKD 
PROTOCOL 

The complete QKD protocol with four-dimensional po- 
larization states exploits five mutually unbiased bases 
with four states in each [27j ■ In terms of biphoton states 
the first three bases consist of product polarization states 
of two photons and last two bases consist of two-photon 
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TABLE II: 



protocol 1 


a (deg.) 


theory 


experiment 


F 








1 


-0.0295 - 0.0306i 
0.0543 - 0.0202i 
-0.0154 - 0.0093i 
0.9972 


0.995 


10 


-0.0015 - 0.0229i 
0.0038 - 0.0050i 
0.2013 + 0.1735i 
0.9638 


0.0584 - 0.1926i 
0.0073 + 0.0170i 
0.1633 + 0.1233i 
0.9577 


0.963 


20 


-0.0021 - 0.0386i 
0.0154 - 0.0162i 
0.3430 + 0.3625i 
0.8654 


0.0015 - 0.0326i 
0.0019 - 0.0660i 
0.3967 + 0.2085i 
0.8909 


0.976 


30 


-0.0015 - 0.0445i 
0.0337 - 0.0225i 
0.3530 + 0.5716i 
0.7383 


-0.0310 - 0.0542i 
0.0166 - 0.0416i 
0.4526 + 0.5386i 
0.7065 


0.991 


40 


-0.0005 - 0.0440i 
0.0511 - 0.0116i 
0.1601 + 0.7466i 
0.6421 


0.0112 - 0.0823i 
0.1174 - 0.0328i 
0.1087 + 0.8363i 
0.5167 


0.970 



Here \H) = |1), \V) 
|1>-|0», \R) 



|0), \D) e 
|l) + i|0», 



i 

\L) 



(|X> H- 10», X> = 
l)-i|0» 



^(|1)-|0», |ij) = -^(|l) + i|0)), \L) = ^( 
indicate horizontal, vertical, +45 linear, -45 linear, right- 
and left-circular polarization modes correspondingly, and 
lower indices numerate the frequency modes of two pho- 
tons. It has been proved Q that it is sufficient enough 
to use only first two or three bases for the efficient QKD. 
Using incomplete set of bases one sacrifices the security 
but enhances the key generation rate. Since the realiza- 
tion of Bell measurements for the last two bases requires 
a big experimental effort on both preparation and mea- 
surement stages of a protocol, we will restrict ourselves 
to first three bases. As we will show experimentally, all 
states from the first three bases can be prepared with 
the help of a single non-linear crystal and local unitary 
transformations. This is fundamental distinction with re- 
spect to biphotons-qutrits where SU(2) transformations 
between states from mutually unbiased bases are prohib- 
ited. We also present a measurement scheme that allows 
to discriminate the states belonging to one basis deter- 
ministically thus allowing its implementation in realiza- 
tion of QKD protocol with polarization ququarts. 



TABLE III: 



protocol 2 


a (deg.) 


theory 


experiment 


F 







-0.0555-0.0204i 










-0.0059+0.005i 


0.996 



1 


-0.0425+0.0052i 
0.9973 




0.8097 


0.8067 




20 


-0.4568-0.3527i 


-0.4847-0.3304i 


0.998 


-0.0103-0.0859i 
-0.0316+0.0529i 


0.023-0.0554i 
-0.0174+0.0413i 



TABLE IV: 



protocol 1 


theory 


experiment 


F 


0.707 


0.7326 







0.0818-0.0963i 


0.941 





0.0003-0.0281i 


0.707 


0.6131+0.2657i 




0.707 


0.6597 







0.2518+0.4692i 


0.934 





0.0897-0.0319i 


-0.707 


-0.6155+0.3261i 





entangled states: 



I. Itfrifa); {H^; \VyH 2 )\ \V{V^U 

II. \D 1 D 2 ); \DxD 2 ); j^i-Da); |DiX>a), 

III. \RiB»); \RiL 2 ); \LxRt)\ \LxL 2 ), 

IV. \R 1 H 2 ) + \L 1 V 2 ); \RxH 2 ) - [L^); (39) 

|LiF 2 ) + |iJiV 2 >; |Lii?2>-|#i^2>, 

V. \H 1 R 2 ) + \V 1 L 2 ); \H X R 2 ) - |ViL 2 >; 
|HiL 2 ) + |yii? 2 ); \H x Li) - \ViBa). 



A. Experimental procedure. 

Experimental setup for generation of ququart states 
belonging first three bases 1|39|1 is shown in Fig. 8. To 
varify the state generated in the set-up either protocol 1 
or protocol 2 might be used. Usually we applied protocol 
1 for checking the ququart state. Let us consider as an 



pump 




LiIOj DPI DP? ZP1 ZP2 



FIG. 8: Setup for preparation and measurement of ququarts 
which can be used in QKD. Wave plates DPI, DP2 oriented 
at 45° degrees with respect to vertical axis serve as dichroic 
retardation plate with variable optical thickness which is con- 
trolled by tilting angle 9. Two zero-order plates ZP1, ZP2 
allow to chose the basis. 

example, the preparation of a state \HiV 2 ) from the ini- 
tial state IV1V2). This transformation can be achieved by 
a dichroic wavcplate that introduces a phase shift of 2ir 
for a vertically polarized photon at 605nm, a phase shift 
of it for the conjugate photon and that is oriented at 45° 
to the vertical direction. Using quarts plates as retar- 
dation material it is easy to calculate that the thickness 
of the waveplate that does this transformation should 
be equal to 3.406mm [2^. Since these waveplates were 
not readily available and the result of transformation 
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is extremely sensitive to the small variations of thick- 
ness, we used the following method to achieve the desired 
thickness. Two quartz plates with geometric thicknesses 
of 3.716mm (DPI) and 0.315mm (DP2) with orthogo- 
nally oriented optical axis were placed consecutively in 
the biphoton beam. The consecutive action of these two 
wavcplates correspond to the action of quartz wavcplatc 
with an effective thickness of 3.401mm. If then one can 
tilt these waveplates towards each other by a finite an- 
gle 9, then the optical thickness of the effective wave- 
plate, formed by DPI and DP2 will be changing, and, at 
a certain value of 9, the desired transformation will be 
achieved. This corresponds to maximal coincidence rate 
when the measurement part (protocol 1) is tuned to select 
a state I-H1V2). Monitoring the coincidences, one can ob- 
tain the value of 9 for which the maximum occurs. Then, 
fixing the tilting angle at this value, one can perform a 
complete quantum state tomography protocol in order to 
verify if the state really coincides with the ideal. In order 
to change the basis from / to 11(111), zero order half- 
quarter) waveplates ZPl (ZP2) oriented at 22.5° (45°) 
were used. This procedure is repeated for generation of 
any of the states. Fig. 9. shows the coincidences and sin- 




angle between quartz plates (cleg.) 



FIG. 9: Dependence of single counts (upper) and coincidences 
(lower) on tilting angle 6 

gle count rates versus the change of tilting angle 9 which 
determines optical thickness of the effective waveplate. If 
the measurement setup is tuned to select the state \H1V2) 
then dependence of coincidence rate on the plate optical 
thicknesses Si is given by formula: 



oc (ai t /j 2 t ai^2) = sin 2 (r5i) cos 2 (t) 2 ), 



(40) 



whereas the single counts distribution in the upper chan- 
nel is given by 



l702nm OC (ai^Ol) = Sm 2 (r5i), 



(41) 



where optical thickness depends on tilting angle as fol- 
lows: 



x nh i 



{n ej f 



K) 2 



sin* 



sin 2 9 



(42) 



The solid lines in Fig. 9 show the theoretical curves. 
We performed tomography measurements for the both 
maxima as well as for the minimum. The minimum in 
coincidences occurs when intensity in any channel drops 
to zero, so it is not a necessary condition for distinguish- 
ing the orthogonal state to that one selected by given 
settings of polarization filters. Nevertheless according to 
calculations and our measurements the minimum in the 
coincidences at Fig. 9 exactly refers to the state \V1H2). 
Starting from the \V702nmV605nm} , we also prepared and 
measured the whole set of states from (I39f) belonging to 
first three bases. The result of state reconstruction is 
presented in Table V. 

TABLE V: 



I 


\H!H 2 ) 


|HiVa) 


\ViH 2 ) 


IV1V2} 


F 


0.98 


0.94 


0.98 


0.98 


II 


\D1D2) 


\DiDa) 


\D X D 2 ) 


|DiDa) 


F 


0.97 


0.96 


0.95 


0.99 


III 


\R1R2) 


\RiL 2 ) 


\l*Ra) 


\I*La) 


F 


0.95 


0.95 


0.96 


0.97 



At the same time the method discussed in this section 
allows one to unambiguously distinguish all states form- 
ing chosen bases. The measurement set-up which has 
been already tested in our experiments is shown in the 
Fig. 10. It consists of the dichroic mirror, separating the 



RNG 



QP HP 



choise of 
basis 



Bob 



DBS 



PBS 



PBS 



H, 



Inputi 



Input2 



Input3 



Input4 



Output 



FIG. 10: Measurement part at Bob's station 

photons with different wavelenghts, and a pair of polar- 
ization beamsplitters, separating photons with orthogo- 
nal polarizations. Four-input double-coincidence scheme 
linked to the outputs of single-photon detectors registers 
the biphotons-ququarts. For example for the first basis, 
the scheme works as follows, provided that Bob's guess 
of basis is correct: 

if state \H1H2) comes, then detectors D4, D2 will fire, 
if state \H1V2) comes, then detectors D4, Dl will fire, 
if state \V1H2) comes, then detectors D3, D2 will fire, 
if state IV1V2) comes, then detectorsD3, Dl will fire. 
Same holds for any of the remaining correctly guessed 
bases, since the quarter- and half waveplates transform 
the polarization to HV basis in which the measurement is 
performed. Registered coincidence count from a certain 
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pair of detectors contributes to corresponding diagonal 
component of the measured density matrix. So if the ba- 
sis is guessed correctly, then the registered coincidence 
count deterministically identifies the input state. We il- 
lustrate this statement by the table which shows total 
number of registered events per 30 sec for the input state 
\R\L2) measured in circular basis and calculated compo- 
nents of experimental (theoretical) density matrix. 



TABLE VI: Coincidence rate and density matrix components 



D 4 D 2 


P11 




p22 


D 3 D 2 


P33 




P44 





0.0(0) 


220 


0.973(1) 


6 


0.027(0) 





0.0(0) 



The main obstacle for practical implementation of free- 
space QKD protocol based on ququarts is that one needs 
to perform fast polarization transformation at the se- 
lected wavelengths. There are different ways how to over- 
come this problem and we will discuss its elsewhere. In 
this section we briefly mention the possible ways. Since 
it is not practical to tune the tilting waveplates every 
time one wants to encode a quart value, we suggest ei- 
ther to use a polarization modulator that operates on 
two wavelengths or to split the photons with a dichroic 
mirrors and perform these transformations on halves of 
a biphoton independently in a Mach-Zchndcr like con- 
figuration. It is important to note that interferometric 
accuracy in Mach-Zchnder is not needed, since it is used 
only for spatial separation of photons. The practical solu- 
tion would be to couple the down converted photons in a 
single mode fiber to ensure a perfect spatial mode overlap 
and then to split them with wavelength division multi- 
plexer (WDM). Then, the switching between the states 
can be done with the polarization modulators that intro- 
duce a 7r or 2-7T phase shifts for the selected wavelength. 
The choice of basis on Alice's side is done by a zero order 
quarter and half waveplates, that can realized within a 



Pockcl (Liquid Crystal) cell driven by randomly selected 
voltage. Free space communication is proposed since it is 
not practical to distribute a polarization state within an 
optical fiber. On Bob's side, the random choice of basis 
(RNG) is performed in a same way as on Alice's side. 
Then photons are spatially separated with the help of 
WDM or dichroic mirror and each of the photon is sent 
to a Brown- Twiss scheme with a polarizing beamsplitter 
that does a projection of an arrived photon on H or V 
state as it shown in Fig. 10. Moreover, registering co- 
incidences allows one to circumvent the problem of the 
detection noise that is common for single-photon based 
protocols. If the coincidence window is quite small, it 
is possible to assure a very low level of accidental coin- 
cidences for the usual dark count rate of single photon 
detectors. 

To conclude we have suggested and tested a novel 
method of preparation, and measurement of subset of 
four-dimensional polarization quantum states. Since for 
this class of states the polarization degree is not invariant 
under SU(2) transformations it is possible to switch be- 
tween orthogonal states using linear transformations like 
geometrical rotations and phase shifts. We discussed the 
family of states that can generate the whole basis states 
using only simple polarization elements and thus might 
be used in useful applications. 

While completing this manuscript we learned about 
closely related work performed by R. B. R.Adamson and 
co-authors |29j . 
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